Electronic correlations on a metallic nanosphere. 
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We consider the correlation functions in a gas of electrons moving within a thin layer on the surface 
of nanosize sphere. A closed form of expressions for the RKKY indirect exchange, superconducting 
Cooper loop and 'density-density' correlation function is obtained. The systematic comparison with 
planar results is made, the effects of spherical geometry are outlined. The quantum coherence of 
electrons leads to the enhancement of all correlations for the points-antipodes on the sphere. This 
effect is lost when the radius of the sphere exceeds the temperature coherence length. 
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Last years a considerable theoretical interest has been 
attracted to the electronic properties of cylindrical and 
spherical nanosize objects. This interest is mostly re- 
lated to the physics of carbon macromoleculesjil where 
the theoretical works are devoted to the band struc- 
ture calculations!! and to the effects of topology for the 
trasport and mechanical properties. The discussion of 
the topological aspects of electronic motion is primarily 
connected with carbon nanotubesu and is based on the 
analysis of the effective models, incorporating the partic- 
ular geometry of the object. 

Other branches of the condensed matter physics, where 
the spherical objects appear, concern the studies af the 
nonlinear optical response in composite materialsu And 
the photonic crystals on the base of synthetic opals In 
case of opals, the Si02 rjanosize balls could be coated 
by metal films as well.El The theoretical studies here 
are focussed on the nonlinear optical properties and the 
metal costing is characterized by an effective dielectric 
functionia This approach should be revised when the 
coating has a width of aijew monolayers ,0 a limit allowed 
by modern technologies a 

Thus the investigation of the electronic properties per- 
tinent to the spherical geometry can provide an interest- 
ing link between the carbon macromolecules and other 
spherical structures. On the theoretical side of the prob- 
lem, one meets a unique possibility to establish a bridge 
between the methods of microscopic many-body the- 
ory and the intrinsic spherical geometry of the atomic 
physics. To be consistent, these theoretical efforts should 
comprise the systematic comparison of the quantities 
found for the planar and the spherical geometries of the 
electronic motion. In a recent paperO we investigated the 
electronic gas on a sphere in a uniform magnetic field. 
The exact solution of the problem was found, two phys- 
ical effects were predicted. First is the jumps in the 
magnetic susceptibility at half-integer numbers of flux 
quanta, piercing the sphere. Second is the localization of 
the electronic states near the poles of the sphere at high 
fields. 

In this paper we study the correlations in the elec- 
tron gas moving on the surface of the sphere. We ob- 
tain a close form of expressions for the indirect RKKY 



exchange, superconducting Cooper loop, density-density 
correlator. The effects peculiar to this geometry are elu- 
cidated, their meaning for subsequent theoretical inves- 
tigations is discussed. It is particularly shown that the 
coherence of the electronic motion on the sphere results in 
the enhancement of the correlations between the points - 
antipodes on the sphere. This effect is lost with increas- 
ing the temperature T, when the temperature coherence 
length becomes less than the radius of the sphere. 

We consider the electrons moving on a surface of the 
sphere of radius tq. The Hamiltonian of the system is 
given by TL — — ^ — h U(r), where m e is the (effective) 
mass of an electron and we have set % = c = 1 . The total 
number of electrons N (with one projection of spin) is 
fixed and defines the value of the chemical potential fi 
and the areal density v = A r /(47rrg). We assume that 
the potential U(r) confines the electrons within the thin 
spherical layer Sr <C ro, i.e. U{r) = at tq < r < ro + Sr 
and U(r) —* oo otherwise. The radial component R(r) of 
the wave function is a solution of the Shrodinger equa- 
tion with the quantum well potential. We adopt that the 
chemical potential \i lies below the first excited level of 
R(r), which means 8r < v~ x l' 2 . It possible then to ignore 
the radial component and put r = ro in the remaining 
angular part of the Hamiltonian 7ir> . The eigenfunctions 
to the Shrodinger equation Tin^ — E^> are the spherical 
harmonics Y/ m and the spectrum is that of a free rotator 
model : 

*(0,0) = r o - 1 r /m (0,0), E l = {2m e r 2 )- 1 l(l + l). (1) 

According to (|l|) ^{6,4>) is normalized as 
7-q J | 'J | 2 sm QdOdcf) = 1, it facilitates the comparison of 
our results with the case of planar geometry. 

For the two points r (6, 0) and r' (6' , cf)) we de- 
fine the distance f2 on the sphere as cos f2 = cos 9 cos 9' + 
sin 9 sin 9' cos cj>. One can findQ an exact representation 
of the electron Green's function through the Legendre 
functiont2l 

777- 

G(Sl,ioj n )=-— ^P_ 1/2+a (-cosfi), (2) 
2 cos ira 

where a = y '2m e r^{\x + iu n ) + 1/4 and Matsubara fre- 
quency u> n = itT{2n +1). 
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In the limit a> 1 and for a sin ft > 1 one hasi in the 
main order of a -1 : 

, me C ° S(a7r = °" = (3) 
V27ra sin ft cos 7ra 

At the same time, if a 3> 1 and a(7r — ft) ~ 1 we get 



G(ft ~ tt,w) 



2 cos 7ra 



Jo(a(7r-Q)), 



(4) 



with the Bessel function Jo (a;). When 17 — > 0, the Green's 
function (|J) diverges logarithmically, as it should. 

It is convenient to introduce here the concep t of the 
angular Fermi momentum L — \J 2Tn e v 2 [i +1/4 «-> kpr^. 
For simplicity of subsequent calculations we consider the 
case of integer L, so that at low temperatures the level 
I = L — 1 is completely filled and I = L is empty. It 
is worth noting that the number of electrons N = L 2 
and the energy level spacing at the Fermi level AE ~ 
Anv/{m e L) ~ 2^,/L. For densities ^ ~ 10 14 cm -2 and 
r = 100 Awe have L ~ 30 and AS - 10 meV. 

The oscillating factors exp ±i(an— aft — 7r/4) in (Q) cor- 
respond to coherence of two waves. One propagates along 
the shortest way between two points and another wave 
goes along the longest way, turning around the sphere. 
One can show that this coherence is destroyed by the fi- 
nite quasiparticlc lifetimes] or by the temperature. The 
latter reason takes place at T > AE, in which case the 
Green's function acquires the form (uj n > 0) : 



(5) 



This equation, after natural substitution Lft — * kpr, cor- 
responds to the usual expression for the planar geometry. 

Let us now discuss the magnetic correlation function 
on the sphere. This quantity defines the indirect RKKY 
exchange interaction between the localized moments and 
thus is related to the transverse NMR relaxation rate. 
It is also reponsible for the magnetic instability in the 
electronic system (spin-density wave state). The range 
function of the RKKY exchange interaction is given by : 



X 



(n) = -Tj2G 2 (n,iuj n ). 



(6) 



We discuss first the case T — » and use the limiting 
relation TJ2 n ~ > /f^, duj/(2n). Using (Q) we write 



L-l 



x(n) 



Aw 2 rn da 



(a+l/2)[P (-cosft)f 



1=0 



(8) 



This expression was first obtained by Larsen0 in a 
different way. Being formally exact, the expression (JsJ) is 
not particularly useful at L ^> 1 . Jn the case of large L 
and at ft ~ it it was suggested inO to approximate the 
sum (0) by an integral, thus obtaining 



-L[P L {- cos ft)]' 



(9) 



This estimate is however unsatisfactory at ft ~ 1 , since 
at large v one has P u (cosQ) ~ v~ x l 2 cos(z/ft — tt/4). 
Then the leading terms in the sum (||) acquire the form 
J2i cos 21Q, and the consideration of the next-order terms 
is needed. 

To find x(f2) in the closed form at L ^> 1, we return 
back to (||) and use the asymptotic expressions (||), (Q) 
for the Green's function. Next we rescale u> — > 2/iuj/L 
(it corresponds to measuring the energy in units of spac- 
ing AE) and neglect the terms oc l/L in the expansion 
a = \/L 2 + 2iLuj ~ L + iuj + uj 2 /2L. After that the inte- 
gration over uj in is easily done and we come to the 
improved estimate for (||) in the form : 



X(«) 



4tt 



47r 3 sin ft 



-n)], 

ft-7T 

sin ft 



(tt-Q) < l/L 



sin(2Lft) 



ft ~ 1, 



(10) 



(11) 



with the Bessel function Jo(x). One can verify that the 
Eqs. @ and (|TTJ) match smoothly at (tt-CL) ~ l/L. It 
is convenient to rewrite ( [Tl|) further as x(^) — X^(^) + 
X (2) (^) where 



x (1) (^) 
x (2) (^) 



4tt 2 



~ 2 sin(2Lft) Q 



sin 

-2 



2 ft 



ft 



47r 3 sin ft 



7T 

2 

sin(2Lft) 
sin ft 



(12) 



ft -tt9 



„ 71" 

n "2 



and the step function Q(x) = I &t x > Q. For the small 
angles l/L <C ft <C 1 the x^ term dorpjaates and we 
recover the result of the planar geometryMEj 



x(fi) 



da- 



cos 7ra 



[P a -i/ 2 (-cosft)] 2 (7) 



where the contour C in the complex plane starts at 
passes through the point L on the real axis 
and goes to e l7T ^ 4 oo. We can shift the integration con- 
tour to the imaginary axis of a, where the integral is zero 
due to the oddness of the integrand in (Q). Upon this 
shift we intersect the double poles at half-integer real a. 
Therefore the above integral reduces to the finite sum of 
residues and is equal to 



(27rr ft) : 



■ sin(2Lft) 



(27rr) : 



■sin(2fc F r) (13) 



Discussing the analogue of the q— representation of 
(|li"|), we note that the usual formula for the inhomo- 
geneous static susceptibility x(°±) = / dre lqr x(r) is re- 
placed by Xlrn = r Q 

J x(^)YffMl, if) sin ft dVtdip. In the 
absence of tp in x(^)j we haveEJ : 



Xi = 2-nrl I dft sin ft P;(cos ft)x(ft) 



(14) 
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Recalling the property Pi(— x) = (—l) l Pi(x) and noting 
that x^H 71 — fl) = X^(^)i we conclude that the term 
X^(&) contributes to all harmonics Xh whereas the "co- 
herent" part x^(ty contributes only to xi with even 
Z. Below we show that x^ disappears at high enough 
temperatures, T ~> AE, or in the presence of scattering. 

Calculating the contribution of the term x : we en- 
counter the integral J^ 2 dD.Pi (cos fl) sin(2LJ7)/ sin(fi) 
which is determined by small fl. Thus one can use 
Pi(cosfl) ~ Jo[(l + l/2)f2] and extend the integration 
to the infinite upper limit. 

Further, the first term in x^ yields the integral 

f* dnp 2 i(cosn) = r 2 (z+i/2)/r 2 (z+i)~ (z+i)- 1 . The 

second part of x^ (fl) proportional to fl sin(2Lfl) / sin 2 fl 
is also important at fl <C 1 ; we can approximate the cor- 
responding integral by L dflJo((l + ^)fl) sin(2Lfl). 
Combining all the contributions, we have 



Z7T 



J_ ( 1 + 1/2 
irL 12 V 2L 



Xl 



— f (— 

2ir J1 \ 2L 



for even and odd Z, respectively. Here the function 

f 1 (x)=7r/2, x<l 

= sin _1 (l/x), x > 1 



(15a) 
(15b) 

(16) 



corresponds to x (^)* K defines a cusp in at 
Z = 2L ^-t—.2kp and thus resembles the usual planar 
expression.!! 2 ] The function 



f 2 (x)=x~ 1 + (l-x 2 )~ 1 / 2 , x<l 



(17) 



x > 1 



stems from the part %( 2 )(r2) and is present in the even 
harmonics x.i- This coherent term has a prefactor 1/L 
in ( |l5| ) and formally vanishes in the limit L oc ro — > oo. 
However, /a (a;) is singular at x = and at x = 1, which 
points deserve more attention. At Z ~ 2L one estimates 
the contribution of 

x (2) „ L -l/ 2j j e 

much less than 

X^ 1 ) ~ 1. At Z ~ 1 both contributions are comparable. 
Moreover, using (|§|) and relevant formulas ircj one can 
find the exact identity for the uniform static susceptibil- 
ity xi=o — 0- The vanishing value of xi=o is connected 
with our choice of the chemical potential lying between 
the discrete energy levels L— 1 and L, so that the excited 
electronic states are separated from the vacuum ones by 
the energy AE. It resembles the situation, with vanishing 
static susceptibility in superconductors.E3 Exploring fur- 
ther this analogy, it is interesting to observe that in the 
clean metal and at T — > the coherent term f 2 provides 
a maximum of xi a t Z = \/2L. This poini— .similarly to 
the case of ferromagnetic superconductoral3'E£l, could be 
labeled as the momentum of helikoidal magnetic order- 
ing. 



Next, we discuss the superconducting correlation func- 
tion, which reflects the tendency to superconducting pair- 
ing and defines on the mean-field level the superconduct- 
ing temperature, T c . The basic object here (found also in 
the localization theory) is the static Cooper loop, which 
could be written in r— representation as 



n(ft) = T$^G(n,»w n )G(n,-tw n ), 



(18) 



Unlike the above case of the magnetic loop xfl)> we 
could not arrive at some closed type of expression for 
n(f2) with the use of the exact formula (||) for G(fl,u>). 
However, one can find an approximate form of II(fi) ex- 
ploring the large L asymptote (|^) for the Green's func- 
tion. A calculation similar to the above one gives the 
result (cf. (|ll])) 



rr(fi) 



Air 2 



■LJ$[L(w- fl)} , (tt - fl) < 1/L (19a) 



47r 3 sin fl 



fl 



sin fl 



- sin(2Lfl) 



fl ~ 1 (19b) 



Again, Eqs. (|19|a) and (|19|b) smoothly match at (tt — 
fl) ~ 1/L. As above, we decompose Yl(fl) into the sum 
II«(fi) + U^(fl) with 



nW(fi) = 

n^ifi) = 



vnr. 



„-2 



e 



47r 2 sin 2 fl 

2 r TrG[fl 



2 



(20) 



mr, 



47r 3 sin fl 



tt/2] - fl 



sin fl 



sm(2Lfl) 



The property n( 2 )(7r - fl) = -Tl^fl) should be noted 
here. At 1/L < fl < 1 the term n (1) (£!) is a principal 
one and the planar resultO is restored 



U(fl) 



An 2 sin fl 



(2ttt) 2 



(21) 



Finding the momentum representation of (|l9|b) accord- 
ing to (pT[), we encounter the logarithmic singularity at 

— ^ and all Z, which originates from the term nW(f2). 
This divergence should be cutoff at the lowest allowable 
angles flo ~ ao/ra with the interatomic spacing ao. 

At the same time, the part II' 2 ' (fl) provides a conver- 
gent integral which contributes only to odd harmonics, 
in view of the abovementioned property. Further, it sat- 
urates at small fl, where H^(fl) coincides with x^fl). 

The whole result is then given by : 



n, 



2tt \aoZ / 



' 2tt 



1 a l J ttL j \ 2L 



(22a) 
(22b) 



for even and odd Z, respectively. Comparing it to ( |l5| ) we 
see that the same coherent term f 2 (l/2L) contributes the 
harmonics of different parity. A simple analysis shows 
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that the presence of fo(lj2L) in (|5h does not modify 
an overall character of 11/ for all L> 1, and it is the 
logarithmic term that defines the tendency to the super- 
conducting pairing at / — > 0. It should be also stressed 
that the static Cooper loop has a singularity at an analog 
of 2kp, the feature absent in the planar geometry. 

Let us discuss now the 'density-density' correlation 
function of the electronic gas on the sphere. This function 
describes the variation in the electronic density, caused 
by impurities, and can be written asu3 

C(r, r') = ~((n(r) - v)(n(r') - v)) - S(r - r') (23) 



= -\ (r^G(r,r^0e-^ 



where r — > +0 in the last line. 

In the planar geometry we have at T - 



(24) 



with T(x) — x/smhx. Therefore the temperatures 
T > AE lead to the exponential decrease of the cor- 
relations at > AE/(irT) = £t/?~o, with the tempera- 
ture coherence length £y = 2/i/ (irkpT). One should also 
note the disappearance of the coherent terms x (ty an d 
II^ 2 ) (Q) present at Q <C 1 in magnetic and superconduct- 
ing correlators, ([l2]) and (|20|), respectively. 

Summarizing, we find a closed form of the correlation 
functions for the electron gas on the sphere. The effects 
peculiar to this geometry are analyzed, the role of finite 
temperatures for our results is elucidated. 

I thank S.L. Ginzburg, S.V. Maleyev, A.G. Yashenkin, 
A.V. Lazuta for stimulating discussions. The partial fi- 
nancial support from the Russian State Program for Sta- 
tistical Physics (Grant VIII-2), grant FTNS 99-1134 and 
grant INTAS 97-1342 is gratefully acknowledged. 







CM 



1 

nr 



sin(2fci?r) — 1 
n 2 kpr 3 



(25) 



with r = |r — r'|. We see that at large r the Friedel 
oscillations with a period (2kp)~ 1 take place. 



For the spherical geometry and large L we use (|3|), (|24|), 
to obtain in the limit T = : 



CM 



1 sin(2Lft) - 1 
47r 2 r 2 Lsinf7sin 2 (fi/2) 



(26) 



This expression holds at L sin fi > 1 and has the obvi- 
ous correspondence with ( |25| ) for fi < 1. At large il ~ ir 
the amplitude of the correlations grows, in accordance 
with the previous formulas (|TT|), (fl9|). 

Now let us discuss the case of finite temperatures. 
First, we note that for the considered case of the fixed 
number of particles, the position of the chemical poten- 
tial depends on T and is determined from the equation 
N = ^2 l (2l + \)np(Ei), with Fermi function np (x). Some 
analysis shows, however, that the value of fj, at finite T 
does, not deviate much from its zero-temperature value 
/ioO Specifically, this deviation satisfies the inequalities 
\fi — fJ@l/T < 1 and \fi — (j,q\ < AE/2 and is ignored 
below.El 

The above correlation functions are not essentially 
modified at T <C AE. The changes occur when T > AE, 
in which case we use Eq. (§) and perform the Matsub- 
ara sums in (0), (|l8|), (pi}). As a result, we arrive at the 
following expressions : 



. msin(2LO) /2vrTO\ 
X(0, T) ~ — 



(27rr fi) 2 \ AE J 



n(n,r) 



-T 



/27rra\ 



47r 2 ^ 2 V AE 



■ 



7T 2 r 2 L0 3 
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If the Fermi level L — 1 is not completely filled, fio is defined 
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as the limiting value of /i at T — > 0. the fact that in the planar limit AE — > we recover the 

The relative smallness of this deviation is connected with two-dimensional, and hence constant, density of states. 
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